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Abstract 



A renormalization group (RG) analysis of the superconductive instability of 



, an anisotropic fermionic system is developed at a finite temperature. The 

^5 method appears a natural generalization of Shankar's approach to interacting 

a; 

, fermions and of Weinberg's discussion about anisotropic superconductors at 

Q I T = 0. The need of such an extension is fully justified by the effectiveness 
of the RG at the critical point. Moreover the relationship between the RG 

^ \ and a mean- field approach is clarified, and a scale-invariant gap equation 



is discussed at a renormalization level in terms of the eigenfunctions of the 
interaction potential, regarded as the kernel of an integral operator on the 
Fermi surface. At the critical point, the gap function is expressed by a single 
eigenfunction and no symmetry mixing is allowed. As an illustration of the 
method we discuss an anisotropic tight-binding model for some classes of high 
Tc cuprate superconductors, exhibiting a layered structure. Some indications 
on the nature of the pairing interaction emerge from a comparison of the 
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model predictions with the experimental data. 
PACS numbers: ll.lO.Hi, 64.60.Ak, 74.20.-z, 74.72.- 



I. INTRODUCTION 

In a recent review,0 Shankar unified a large variety of pliysical plienomena and tlieir 
conventional descriptions, such as Landau theory of Fermi liquids, charge density waves 
and nesting in the Fermi surface, superconductive instabilities and mean-field theories etc., 
under the scheme of a renormalization group (RG) approach, with the aim of a better 
understanding of the general nature of the physics involved in such processes. Central to 
that approach is the idea that the instability towards a phase transition in an interacting 
fermionic system may be discussed in a very general fashion, using the RG technique, which 
involves the cut-off and the rescaling of the energies and the momenta close to the Fermi 
surface. For instance, a charge density wave or a superconducting condensate are recovered 
owing to very general symmetry properties of the single particle dispersion relation, i.e., the 
existence of a nesting vector and the validity of time reversal symmetry, respectively. 

In an almost contemporary workjl Weinberg derived a RG flow equation for supercon- 
ductors whose Fermi surface satisfies just time reversal invariance, employing the standard 
field-theoretic scheme. In particular, Weinberg's work doesn't claim for a spherical Fermi 
surface, and therefore shows itself more suitable for taking into account anisotropic materials 
such as the high layered cuprates. 

An extension to a non-zero temperature of the RG analysis is called for, at least for 
two major reasons: i) the RG approach is more effective around the critical point, which is 
generally located at a finite temperature T^, ii) in order to make contact with the conven- 
tional mean-field approach and with the experimental data, the RG predictions would be 
required at a finite temperature too. Such an extension is straightforward, from a technical 
point of view, and even for a generic anisotropic system, the flow equations are decoupled 
in terms of the eigenvalues of an integral operator V, whose kernel is the marginal pairing 
coupling evaluated at the Fermi surface.i The critical temperature Tc is then found out to 
be a simple function of the most negative eigenvalue, which determines the transition point. 
The standard BCS expression is recovered for a spherical Fermi surface with a constant 
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pairing coupling, whereas for a generic rotationally invariant system the eigenfunctions re- 
duce to spherical waves, since then the angular momentum and V are mutually commuting 
operators. 

In order to make contact with the standard mean-field approximation, we recover again 
the same analytical expression for T^, though following a quite different path, starting from 
the usual gap equation. The occurrence in such an expression of just one eigenvalue is a token 
of the existence of a leading part in the pairing coupling, which leads the transition. In fact, 
the expansion of the integral operator V in terms of its eigenfunctions allows one to determine 
its relevant part as its projection over the eigenfunction belonging to the most negative 
eigenvalue. In other words, close to the critical point, the interaction operator V acts on the 
one-dimensional subspace generated by just one eigenfunction. All the physical quantities 
depend on this eigenfunction, which entirely determines even the symmetry pattern of the 
gap function. As a consequence, the occurrence of any symmetry mixing in the gap function 
(such as s-d wave mixing) is ruled out at the critical point, since in that limit the gap 
function reduces to a single eigenfunction with a fixed symmetry (provided that the most 
negative eigenvalue is not degenerate). 

Even far from the critical point, the expansion of V in terms of its eigenfunctions may 
be inserted in the gap equation, thus providing a scale invariant relation, which is of some 
utility especially when dealing with separable potentials, when the number of eigenvalues is 
finite and so is the expansion in terms of the eigenfunctions. The resulting set of coupled 
non-linear equations always admits symmetric solutions, i.e. solutions for the gap function 
which share the eventual symmetry of the physical system. However, far from the critical 
point the equations are highly non-linear, so that the uniqueness of the symmetric solutions 
is not guaranteed, and broken-symmetry solutions are not forbidden. 

Such general aspects of the solutions, as their symmetries near the critical point, are 
more easily discussed by an RG approach than within a mean-field approximation. This 
feature, together with a major handiness in dealing with the numerical cases, provides the 
RG approach with more appeal than the mean-field approximation, without spoiling their 
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equivalence. 

A nice illustration of the method is provided by a simple tight-binding model, recently 
proposed by Spathis et a/.i^ for the single-particle dispersion relation of a layered high 
superconductor, namely Bi2Sr2CaCu208+5 (BSCCO). In order to compare our results with 
those already obtained with a conventional BCS procedure,!'^ the same dispersion relation is 
exactly employed, where rotationally non-invariance accounts for the structural anisotropy in 
BSCCO, and two singlet pairing couplings are considered, namely an on-site and a nearest- 
neighbour interaction. 

In particular, the characteristic dependence of Tc versus the carrier concentration N is 
recovered,^ although it seems to be a mere effect of the peaked quasi-bidimensional density 
of states. Nonetheless, we observe the best qualitative agreement with the experimental 
data when the two couplings have opposite signs, and precisely in the case of an attractive 
on-site interaction (negative HubbardI U) and an inter-site nearest-neighbour repulsion. 

Moreover, while a nearest-neighbour attraction always yields a superconducting ground 
state at T = 0, for any value of the chemical potential fi (though the critical temperature 
may be extremely low), an on-site attraction is cancelled out by the presence of a nearest- 
neighbour repulsion, for a quite large range of values for the chemical potential. 

From a physical point of view, a negative Hubbard U could be the effect of a very short 
ranged coupling interaction, while a nearest-neighbour repulsion could be justified by the 
Coulomb long-range potential. 

We observe, however, that such comments on the nature of the pairing interaction are 
far from being conclusive, and a larger number of pairing couplings should be retained in 
the model potential for a full comparison with the experimental data. From this point of 
view, the method is easily implemented, since adding any other coupling merely implies an 
increase in the space dimension of the integral operator V. 

The paper is organized as follows. In Section II we outline the RG approach at a finite 
temperature for a generic anisotropic fermionic system. In Section III very similar results are 
recovered by use of a standard gap equation, and the possible predictions of the RG approach 
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about the structure and symmetry of the gap function are discussed. Later in Section IV, 
as an illustration of the method, we consider a tight-binding model for the single-particle 
dispersion relation in layered high Tc superconductors, and finally in Section V the numerical 
results are discussed. 



Shankar'su RG approach to the superconductive instability properties of an interacting 
fermionic system is here generalized to a finite temperature and for a generic anisotropic 
Fermi surface. No other special symmetry is assumed than time-reversal. 

The RG flow for anisotropic superconductors has been first discussed by Weinberg! at 
zero temperature, whereas the RG analysis turns out to be more effective around the critical 
temperature T^. An extension to finite temperature of Shankar'^ derivation of the flow 
equations is quite straightforward, and we shall focus on the main aspects, thus referring to 
Shankar's pedagogical revie for the details. 

The partition function Z of a many-fermion system can be expressed as the functional 
integral: 



being k a reciprocal lattice vector, uJn = (2n + l)vr//5 a fermionic Matsubara frequency, 
(3 = the inverse temperature (hereafter, we set % = = 1). The part of the action 
which is quadratic in the fields accounts for the free evolution of the system, and may be 
written as: 



II. RG FLOW AT A FINITE TEMPERATURE 




(1) 



being S = Sq + S\rA the action as a functional of the spinor Grassmann fields 
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being £(k) the free single-particle dispersion relation. 

The most general two-particle interaction in d — 3 dimensions is seen to contribute to 
the action through the term: 

-^int = ^ / d(4)d(3)d(2)d(l)5(4321)^(4)^(3)«(4321)^(2)^(l), (4) 

being 

/dW^^EE/^, m^i^aAKu^n,) (5) 

and 

^(4321) = (27r)=^/35(3)(k4 + kg - k2 - ki)5„,+„3,„,+„„ (6) 

where the Dirac S^^^ enforces momentum conservation up to a vector in the reciprocal lattice. 
The interaction m(4321) = Ua4^a3a2ai0^4, ks, k2, ki) satisfies the general symmetry properties: 

m(4321) = -m(3421) = -m(4312) = m(3412) (7) 

and acts asa2x2x2x2 array in the spinor space. 

Central to the RG approach is the assumption that all the integrals may be evaluated 
within the shell |£(k) — n] < A, with A a proper energy cut-off, provided all the parameters 
in the model are renormalized, in order to take into account the elimination of the modes 
with |£(k) — nl > A. Thus we suppose that all such parameters are renormalized to the 
scale A from the beginning. Of course, a "physical" choice of the cut-off A would allow for a 
comparison of the parameters with the "real world" corresponding quantities. A cut-off in 
energy is preferred, instead of a direct cut-off in momenta, when dealing with a rotationally 
non-invariant Fermi surface. 

When a low-temperature fermionic system is considered, the energy cut-off A can be 
chosen small enough to allow a unique decomposition for any vector k belonging to the tiny 
slice around the Fermi surface |£(k) — //| < A as: 

k = ko + Sk, (8) 
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being kg the vector belonging to the Fermi surface (^(ko) — /i = 0) nearest to k and (5k a 
vector orthogonal to the Fermi surface. Expanding e(k) around kg, only the linear term 
shows to be non-irrelevant in the RG flow, and one can safely write: 

e{k) ^ fi + VF{ko)6k, (9) 

being fF(ko) = |Vk£^|k=ko the Fermi velocity. The cut-off condition is then equivalently 
written as: 

\vF{ko)6k\ < A. (10) 

We may always choose a bidimensional 'vector' of parameters 6 = {61,62) on the Fermi 
surface to individuate kg, together with the energy displacement from the Fermi surface 
e = e(k) — /i = UF(ko)5fc to fix 6k. Such a change of variables is taken into account within 
the integrations according to the rule: 




where 

Vk£-(Vk^iX Vk^2)| (12) 

is the inverse Jacobian function. In practice, we shall omit the e dependence in J, since 
it reveals itself irrelevant in an RG sense, by setting J (6) = J{6,e = 0), which is the only 
marginal term in the expansion of J in powers of e. 
The free action Eq. (^ now reads: 

1 f cP6 d-E 
^o = l^H TTT^J^^) —i)a{6,e,uJn)[iuJn-£\iJa{6,e,uJn). (13) 

An elementary step in the RG flow is next deflned as that transformation on Sq which: 
i) integrates all the modes having |£:| < A/s, with s ^ 1; ii) rescales energies and momenta 
as sojn ^ uj'n, se e', in order to restore the original cut-off; Hi) rescales the flelds as 
5-3/2^ At T = 0, such a transformation leaves the free action invariant, so that the 

latter can be considered as a fixed point. 
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d{6,e) 



d{k) 



At a finite temperature, rescaling the energies as sujn ^ uj'n implies a rescaling of tlie 
inverse temperature itself as jS/s ^ jS'. This is what we must pay for having restored the 
cut-off to its original value. Besides, if any energy scale is associated with the system, such 
as a non-zero temperature or an energy gap, this one undergoes a renormalization flow, 
while the cut-off A is kept fixed. If one identifies s = e*, then: 

^ = -A ^ Pt = f3oe-\ (14) 

which shows that renormalizing the action at a fixed energy cut-off A yields a larger effective 
temperature. 

A completely equivalent approach consists in keeping fixed all the energy scales (tem- 
perature, gap etc), while assuming a flow in the cut-off as: 

At = Aoe-K (15) 

One approach is recovered from the other by a bare change of variables in the integration 
over e. 

Let us now examine how such a RG transformation affects the interaction part of the 
action functional, Eq. (^. At a tree level (Fig. ||), the only marginal couplings which 
survive the RG flow are the lowest order terms in the expansion of m(4321) in powers of Si. 
Moreover, from momentum conservation and phase space arguments, they restrict to only 
two contributions: 

F{ei, 62)= m(4321) with Si = 0, ^4 ■ = ^^2 ■ ^^i, {F coupling), (16) 

V{ei, ^4)= m(4321) with Si = 0, ^4 + ^^3 = ^2 + ^^i = 0, {V coupling), (17) 

where VLi denotes the unit vector in the direction of kj = {6.1,81 = 0). 

We are interested in the flow of the \^-couplings which are responsible for the supercon- 
ductive instability of the Fermi surface. We assume the interaction to be spin independent, 
and we take in consideration only spin singlet pairs, so that we can get rid of all the spin 
indices, by posing: 
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V{ei, 64) = ¥{64, 0i)= ^[Mtut(k4, -k4, -ki, ki) - Mtin(k4, -k4, -ki, ki)] = 

= ^['«iTTi(k4, -k4, -ki, ki) - M|-ri|(k4, -k4, -ki, ki)], (18) 

being the u interactions restricted on the Fermi surface. Besides, the general symmetry 
properties, Eq. (^, imply that V{9i, 9^) is symmetric with respect to the inversion tli —tli 
or fi4 — > — ^4. 

Following analogous arguments as for T = ojll we observe that at one-loop level only the 
BCS diagram contributes to the flow in V (Fig. |^) with: 

dr (^1, ^4) = -\l d(6)d(5)G°(6)G0(5) V(^^i, e^)V{e^, ^4)5(6521), (19) 

Z JdA 

where G^{i) denotes the free fermion propagator, which at a finite temperature may be 
written as: 

G'{^) ^ G'iK co^J = , ^^-T^, (20) 

and where /^y^^ d(6)d(5) does not contain the sum over spin indices, and the integration is 
performed over the tiny slice |e(k) ± A| < |dA|. 

For a generic anisotropic Fermi surface, when an interaction is switched on, both the 
chemical potential and the shape of the Fermi surface change according to the definition:!^ 

/i'-£(k) -ReS(k,/i') = 0, (21) 

being /i' the chemical potential in presence of the interaction and S(k, cj) the self-energy 
in the interacting fermion propagator. At one-loop level, a contribution to S comes from 
the 'tadpole' diagram shown in Fig. |l|. However, the propagator in the one-loop correction 
Eq. (|l^) is a free propagator, and its denominator in Eq. (pO|) is determined in terms of the 
free single-particle dispersion relation and chemical potential. So, no serious problem arises 
at one-loop from the renormalization of the Fermi surface. 

In Eq. ([T9|), the incoming momentum and frequency are zero, and therefore the loop 
momenta and frequencies are opposite. From time reversal, e(k) = £:(— k), and taking into 
account the 6 constraint, one gets: 
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where we made use of the Mittag-Leffler expansioniB 

iV-^^ = -tanh^. (23) 
Integrating at the cut-off, with dt = |dA|/A, one finally obtains: 

^-^-(^)/(&^'^'^'<*"^'^''*'^^^' 

where Pt flows according to Eq. (|l^). 

We notice that the choice of a cut-off in energy allows one to take the factor tanh(A/?t/2) 
out of the integration over the Fermi surface, like in the rotationally invariant case. For 
/3o — ^ oo, Eq. ( p4|) reduces to the flow equation for an anisotropic superconductor recovered 
by Weinberg! at T = 0. Besides, T increases monotonically as t flows down to oo, so that 
T — > oo (/5 = 0) is always a fixed point, since Eq. p^ ) would then yield dVt/dt = 0. 

Introducing the measure 

one may regard Eq. ( Jl^ j as the flow equation for the integral operator: 

(9) = I dTe,Vti9,9'M9'), (26) 

where the coupling Vt{9, 9') plays the role of a symmetric Hilbert-Schmidt kernehi Eq. ( plj) 
then reads: 

The integral operator Vt is Hermitean and admits a complete set of orthonormal eigenfunc- 
tions $Q,: 
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(Vt^a) (9)= \tia)^a{0) (28) 
J dr,<l>:(0)<l>^(0)=(5,^. (29) 

As suggested by Weinberg,! Eq. ( p?]) implies that [dV^, Vt] = 0, so that the eigenfunctions 
$Q,(^) do not flow, whereas the eigenvalues do flow according to the (decoupled) equations: 

^M!)=_t,„h(M) (30) 

which can be integrated together with Eq. ([T^ to yield: 

^* ~ 1 + Ao /o tanh (^e-) dr ' ^^^^ 

being Aq the unrenormalized eigenvalue. The fixed point at /? = is generally reached unless 
Xt diverges for some t, where the flow stops. This may only occur if the unrenormalized 
eigenvalue were negative at the beginning of the flow, Aq < 0. As shown in Fig. ^ at variance 
with the T = analysis, this instability may not occur if Aq is small, and a suitable definition 
of the critical temperature can therefore be obtained from Eq. (^) in the limit i —>■ oo, 
i.e., 



oo 



1 + Ao / tanh i^e-A dr = 0, (32) 
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or, after a change of variables. 



1 f^^c ^ 

— = -/ —tanh-, (33) 
Ao JO 4 ^ 

with Ao the negative eigenvalues being largest in modulus, i.e. the first leading to a diver- 
gence in Aj. 

Eq. ( |33D is quite familiar, and in the limit Af3c 1 niay be analytically approximated 

as: 

T,= ^Aei/^«, (34) 
being 7 ~ 0.5772 the Catalan-Euler's constant. 
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We recover the BCS limit for a constant coupling V = —2v, since then: 

Ao = —2v J dTQ = —nv, (35) 

being n the density of states per spin at the Fermi energy. In this limit, a natural cut-off 
is provided by the Debye frequency. Moreover, Eq. is very general, and holds for any 
anisotropic superconductor. The critical temperature is determined by the most negative 
eigenvalue of the coupling kernel Vt{6,6'), while the remaining eigenvalues do not play any 
role in the proximity of the critical point. This idea will be developed in the next section, 
where the role of the eigenf unctions will be clarified in connection with the scale invariant 
gap equation. 



III. MEAN FIELD VS RG: A SCALE INVARIANT GAP EQUATION 

In the mean-field theory of superconductivity, one is usually led to consider the energy 
gap function A(k). It is a quantity of central interest, since it is directly measurable and its 
symmetry patterns are strongly related to the nature of the pairing, as we shall see in the 
following. 

It would be desirable that the RG analysis could give some indication about the nature 
and the size of the gap function. On the other hand, we would like to throw some light on 
the connection between a RG approach and a conventional mean-field approximation. 

At the end of the previous Section, we noticed that only the most negative eigenvalue 
of the integral operator Vt is responsible for the location of the instability. We are therefore 
led to expect that the associated eigenfunction should determine the gap structure at the 
transition point. Indeed, expanding the integral operator Vt in terms of its orthonormalized 
eigefunctions: 

{Vt ■ ) {e) =J2^:{e)Xt{a) [ dre^^^e') ■ , (36) 
we may decompose the coupling kernel as Vt{9, 0') = V^{0, 0') + Vt{6, 9'), where: 
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(37) 



being At(0) the most negative eigenvahie of Vf The divergence of Af(0) at the transition 
point as t — s> oo clearly indicates the relative importance of Vf in driving the instability in 
the proximity of the transition point. In other words, at the transition point, the kernel 
Vt{9,9') may be approximated with its 'critical' part Vf{9,9'), whose structure is entirely 
determined by the eigenfunction $o(^) belonging to Xt{0). 

The very same conclusion may be reached through a completely different path, starting 
from the standard mean-field gap equation:li 

1 



A k = — / — — -V(k,k')—-^t&nh'—^^ 
^ ' 2 J 2 27r 3 ^ ' ^^ k' 2 



(38) 



where V{k, k') is understood as an extension of the kernel V{9, 9') out of the Fermi surface. 
If A is small, compared to the other proper energies of the system, we may assume that 
the important contribution to the integral Eq. (|38|) comes from the region around the Fermi 
surface. Therefore, introducing a cut-off A ^ A, substituting the renormalized Vt kernel for 
the unrenormalized interaction, and neglecting the energy dependence of A and V, we may 
eventually write Eq. (R3) as: 



Md) = -\j dTe.Vt{9,9')At{9') f^de 



tanh 



/3 



(39) 



An explicit t dependence had to be attached to A and (3, albeit trivial, being a conse- 
quence of choosing to keep A fixed during the flow: both A and are energy scales of 
the system, and are therefore expanded as A^ = Ae*, Pf'^ = /3~^e* during the flow. In the 
present context, it would be preferable to work with a flowing cut-off A^ = Aqc^* and fixed 
values of the energy scales A, (3, so that changing variables e i— >■ ee* we may write Eq. (|39[) 



as: 



A(^) = -- / dTg>Vt(9,9')A(9') f^' de , =tanh 



P 



2 



A2( 



(40) 
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We first discuss Eq. (1^) near the transition point, P Pc, where the RG approach 
shows itself more effective than at /3 ^ Pc- In a mean-field framework, the critical point 
would then be defined just by the vanishing of the gap function. Without loss of generality, 
we may pose: 

A(^) = A,xW, (41) 

being x(^) regular at the critical point and Ac a scale parameter, Ac = 0{13 — j3c). Insertion 
of Eq. (^iD into Eq. (^Of) at the critical point yields: 




which shows that x(^) is the eigenfunction of the integral operator Vt belonging to the 
eigenvalue At such that: 

^ = - ^tanhi. 43) 

The latter result makes sense only if Aj is the most negative eigenvalue of Vt, in which case 
Eq. (^31) is identical to Eq. (p^), thus yielding the same characterization for the critical 
temperature, Eq. (^^. 

As we expected, the functional form of the gap function near the critical point is fixed 

as: 

A{9) = A,%{9). (44) 

We notice that substituting the critical part V^'^ for Vt in the gap equation ( ^OD yields the 
same result for both the critical temperature and the gap structure. The latter finding is of 
particular interest when the eigenfunctions display different symmetries, since only one of 
them is seen to rule over the symmetry pattern of the energy gap, so that no wave mixing, 
such as s-d, may occur as T ^ Tc. 

We remark that Eq. (^) stems from a gap equation linearization near the critical point, 
so that its agreement with Eq. ( ^3|) points forward a substantial equivalence of the RG 
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apporach with the mean-field approximation. However, the former approach displays the 
remarkable advantage of characterizing the transition in terms of the properties of the system 
near the Fermi surface. Namely, the shape of the latter together with the interaction kernel 
evaluated at the Fermi energy entirely determine the RG fiow. In this sense, the RG approach 
allows to classify within a unified scheme the instabilities of fermionic systems sharing the 
same behavious near the Fermi surface, although different in nature. 

On the other hand, integrating modes far from the Fermi surface provides a simplified 
description of the system which still contains all its relevant features. Besides, this descrip- 
tion will prove itself handier when dealing with the numerical implementations, as in the 
next Section. 

In general, inserting the expansion Eq. (^) in the gap equation (^OD, we may formally 
write: 

^(^) = ~J2KiO)\tia) I dTe>^a{e')/\{e') t de , ^ tanh 



(45) 

which means that, at any T < Tc, the gap function is a linear combination of eigenf unctions 
of Vt- This is not a trivial statement, since the eigenfunctions appearing in the expansion 
belong to non-zero eigenvalues, whose number may be finite in several physically relevant 
cases, i.e. when the kernel is separable. 

Eq. (|45| ) is highly non-linear, and may possess more than one solution. For instance, 
broken-symmetry solutions may occur as linear combination of eigenfunctions belonging to 
different invariant subspaces. However, in the limit (3 f3c, the physical solution tends to 
a specific eigenfunction with a fixed symmetry, thus preventing the occurrence of a broken- 
symmetry gap function at the critical point. 

Far from the transition point, the present analysis does not add too much to the general 
comprehension of the problem. For instance, at zero temperature the internal integral can 
be evaluated in Eq. ( ^5]) thus yielding, in the limit A ^ A: 
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2A, 



(46) 



The latter gap equation, valid in the limit A <^ A, can be easily proven to be scale 
invariant according to the flow equation (PUD, and is completely equivalent to the scale 
invariant gap equation recovered by Weinberg.! Moreover, we write the expansion Eq. 
as: 



whose insertion into Eq. 
yields: 



Ai9) = Y.Kid)^H, (47) 

a 

making use of the linear independence of the eigenf unctions, 



A(a) = -Ai(a)5]A(a')(a'|log 



2At 



\a), 



(48) 



being {a'\f{9)\a) = J dre$*,(6')/(6')$Q,(6'), which is a set of non-linear coupled equations for 
the coefficients A{a). Multiplying times A*{a)/Xt{a), summing over a, and introducing the 
scale Aq defined as:i 

■ Ao 



drJA 



iog 



A(^^) 



0, 



(49) 



one eventually obtains: 



Ao = 2Ai exp 



'EJA(«)|VA*(«) 



(50) 



EjA(a)P y 

which is valid only in the limit Aq ^ A^. We observe that far from the critical point 
all the components of the coupling kernel Vt contribute to the determination of the gap 
function. In practice, the present approach has the advantage of dealing with a set of 
coefficients A(q;) whose number may be finite, in the case of a separable kernel. However, 
we cannot escape from a full solution of the set of the non-linear coupled equations (0). The 
present decomposition in eigenfunctions is the generalization of the standard decomposition 
in eigenfunctions of the angular momentum, which are simultaneous eigenfunctions of Vt 
only for a fully rotationally invariant system. 
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If some symmetry invariance is present in both the interaction and in the Fermi surface, 
then the set of equations ( ^S]) may be partially decoupled: in other words, if {$q,(6')} is 
an invariant subspace, and Ug is a unitary representation of the symmetry group {g}, i.e. 
[Ug,Vt] = and {IJ{g)^aid)} = {^a{0)}, Wg E {g}, then we may require that the gap 
function shares the same symmetry: 

U{g)A{e) = e^^(^)A(0), G {g} (|A(0)|2 invariant), (51) 

which is equivalent to saying that: 

U{g)<^a{9) = e^^(3)<l>,(e), V<l>, e {$ J (52) 

(the eigenf unctions belong to the same eig envalue e^^(^)) and A{9) = J^^ A{a)(l>*^{6), 

with $a belonging to the invariant subset 

The set of equations ( ^8]) does admit such solutions with a given symmetry (eigenfunc- 
tions of U) since the matrix elements log |2Ai/A||Q;) vanish if the functions 
belong to different invariant subspaces, provided that |A| is invariant. 

Of course the symmetry invariance of the integral operator Vf does not prevent from the 
occurrence of broken-symmetry solutions, since the non-linearity of the set of equations (^Sj) 
does not guarantee the uniqueness of the symmetric solution. 

IV. A TIGHT-BINDING MODEL FOR LAYERED SUPERCONDUCTORS 

The anisotropic flow equations, derived at a finite temperature in the previous Sec- 
tions, are here employed in the framework of a simple tight-binding model recently pro- 
posed for describing the band structure and superconductive properties of Bi2Sr2CaCu208+5 
(BSCCO).i'@ Alike the majority of the high cuprate superconductors, BSCCO is charac- 
terized by a layered structure, which gives rise to many anisotropic physical properties, and 
further results in a rotationally non-invariant hole dispersion relation £:(k). 

Allowing for nearest-neighbour inter-plane hopping and for nearest-neighbour and next- 
nearest-neighbour intra-plane hopping, in the usual tight-binding approximation, e(k) reads 
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e(k) = A [— 2(cos + cos ky) + 4_B cos k^ cos ky — 2C cos k^] , (53) 

where the components kx,ky, kz of the wave- vector k are measured in units of the respective 
inverse lattice spacings. The constants A,B,C have been determined by comparison with 
photoemission dataBi to be A = 0.05 eV, B = 0.45, C = 0.1 in order to reproduce the 
observed hole density of states for BSCCO. In particular, the condition C <^ B <^ 1 denotes 
a smaller intra-plane next-nearest-neighbour hopping probability than a nearest-neighbour 
one, and an even smaller inter-plane nearest-neighbour one, due to the large plane separation. 
In the following, we shall measure all the energies in units of A. 

In order to set up the model, we must assume the existence of a pairing interaction, 
and the short coherence length would suggest the relevance of the short range part of the 
interaction. Thus, on very general grounds, neglecting the long range contribution, we may 
expand any pairing interaction as the sum of an on-site term and of contributions arising 
from nearest neighbour sites, next-nearest neighbours and so on. 

We do not address any question concerning the physical origin of the pairing interaction, 
but once the free fermion dispersion relation has been fixed on a phenomenological basis, we 
rather wish to explore the main physical consequences arising from the choice of a pairing 
interaction. In other words, we keep the dispersion relation fixed, and change the interaction 
in order to observe the effects on the superconductive instability. 

As an illustration of the method developed before, we here retain only two terms in 
the expansion of the interaction, namely the on-site and the in-plane nearest-neighbour 
singlet pairing couplings. The effects of only an in-plane nearest-neighbour pairing have 
been already considered, together with an identical dispersion relation, in the framework of 
mean-field approximation, by Spathis et a/.@ 

The Fourier expansion of the pairing interaction l^(k, k') reads as: 

V^(k, k') =uo + 2u.j, cos{kx - k'^) + 2uy cos{ky - k' ), (54) 
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being uo the on-site Hubbard U and Ux, Uy the nearest-neighbour interactions in the x 
and y directions, respectively. Enforcing the general symmetry properties Eq. (0) on the 
interaction matrix, the singlet pairing interaction results in the symmetrized part of Eq. ( |5^ 
with respect to the inversion k' — k' or k ^ — k, 

y(k, k')= Uq + Ux[co?,{kx - k'^) + cos{kx + k'^)] + Uy[co?,{ky - k'y) + cos{ky + ky)] = 

= Uq + 2ux cos kx cos k'^ + 2uy cos ky cos k'y. (55) 

The Fermi surface, defined as the locus of the points k in the momentum space verifying 
the equation £:(k) = /x, is a quite complicated, definitely rotationally non-invariant surface, 
whose shape varies furthermore with /x. For = —4 + AB — 2C < ^ < —4 + 45 + 2C = ^2 
it is a closed surface, which for small values of the chemical potential ^ can be approximated 
by an ellipsoid of square semiaxes k^ = ky = {fi + ^i)/{2B — 1), k1 = —{fi + fii)/C, which 
reduces to a sphere since C = 1 — 2B. The Fermi surface becomes an open surface as n 
increases up to the full bandwidth, /X3 = 4 + 45 + 2C, where it reduces to the eight zone 
corners. Fig. |^ shows the Fermi surface within the positive octant of the first Brillouin zone 
for various values of fi. 

It seems convenient to choose 9 = {kx,ky) = {Ox,Oy) as suitable coordinates upon the 
Fermi surface, for each value of fi. Due to its varying shape at increasing fi, coordinates 
6 will be affected by limitations, depending on /i, which we shall keep understood in the 
following, when writing integrals over the Fermi surface. The Jacobian function Eq. (|12D 
may then be straightforwardly worked out as: 

J-\9,e) = ^JaC^ - {fx + e + 2{cos Ox + cos Oy) -AB cos Ox cos Oyf. (56) 

The hole dispersion relation Eq. (|53|) also fixes the density of states of the system, in its 
normal phase, which may be expressed as: 

A further integration from the band bottom fi = yields the total fraction of occupied 
states: 
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de n{e), 



(58) 



together with the 'normahzation' condition at the top of the band, N{fi3) = 1. 

The quantities n(/i) and N{fi) have been numerically evaluated, as functions of the 
chemical potential /i, — /ii < < fi3 (Fig. |§, In particular, n{fi) displays a pronounced, 
yet finite, maximum for ^ = ^2 = —2.0. It may be regarded as a token of the system's 
quasi-bidimensionality, due to its layered structure, and would have been a true van Hove 
singularity, resulting in an infinite peak, if no inter-plane hopping had been considered, 
however small (C = 0). 



For our purposes, the pairing potential Eq. (55) is recognized as the symmetrical kernel 
V{6, 6') in Eq. (0), provided that k^, ky ^ 9^, 9y, k'^, k'y h-> 6'^, 9y, and that 9, 9' are suitably 
restricted on the Fermi surface. This kernel may then be put into the 'separate' form: 



with: 



i,j=0,x,y 



^ uo ^ 
2u^ 
2m„ 



(59) 



(60) 



and: 



u(e) 



' 1 ^ 

COS 9x 
y cos 9y J 



(61) 



Insertion of Eq. (BW) into Eq. (E3) for the eigenvalue problem straightforwardly yields: 



i=0,x,y 



(62) 



Further substitution yields: 



E U^{9) 

i=0,x,y 



A(«)7,(a)- E V^M(^) J ^reUj{9)Uk{9) 



j,k=0,x,y 



(63) 
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which, due to the hnear independence of the functions f/j, yields: 



7fc(a) A(a)(5, 



'ik 



E V^J{UjUk) =0 



(64) 



k=0,x,y 



with: 



(65) 



being, in particular, 2{Ul){n) = n{n) (Fig. ||). 

The linear homogeneous system Eq. (|64D allows then for non trivial solutions '^i{a) {i = 
0, X, y) if and only if the secular condition is fulfilled: 



which allows one to determine X{a) as a function of the parameters uq, Ux, Uy and of the 
chemical potential /x. These should be regarded as the eigenvalues of the operator Vt at 
the beginning of the RG flow (t = 0). Their sign therefore accounts for an instability: the 
existence of a negative eigenvalue A(a) heralds a later divergence pattern for t —>■ tc, tc being 
determined by the most negative of the X{a). The eigenfunctions may then be determined 
by standard algebra through the coefficients 7i(a) of their expansion in terms of the Ui. 

We remark that the main numerical task is to evaluate, once for all, the integrals defined 
in Eq. (|65|) . The secular condition, Eq. (0), can then be easily discussed by changing at will 
the values of the coupling parameters. Besides, the overall formalism allows a straightforward 
generalization, by adding other terms in the development of the potential function, Eq. (^5|), 
thus increasing the order of the matrix (UiUj). 



The results we have exposed thus far may be discussed with respect to the various possible 
values in the couplings, Uq, u^, Uy. 

In the case we may neglect an on-site pairing {uq = 0), Eq. ( |5B[ ) factorizes to yield: 



det X{a)6ik 



E V^AU,Uk) =0, 



(66) 



V. DISCUSSION AND FINAL REMARKS 
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A(0)= 

A(l, 2)= {u. + Uy){Ul) ± ^{u.,-UyY{UlY + Au.Uy{U.UyY. (67) 

Furthermore, in the symmetric case, = Uy, the non-vanishing eigenvalues reduce to: 

X{l,2) = 2u4{U^,)±{U,Uy)], (68) 

which are displayed in Fig. ^ as functions of the chemical potential /i. 

Since (U^) > {U^Uy), for any value of fi, the eigenvalues Eq. (^) are both negative 
for any attractive inter-site coupling Ux < 0, however weak. A superconducting ground 
state is thus predicted for any filling of the band, although the critical temperature is 
undistinguishable from zero even as ^ 0.5. This may be regarded as a mere consequence 
of the very pronounced peak in the density of states (Fig. which determines the behaviour 
of all the eigenvalues. 



As shown in Fig. |^, insertion of the most negative eigenvalue into Eq. (34) yields a 
maximum for Tc, and the existence of a superconducting phase for A^ ^ 0.5, with Ux = —1.0. 
Here, the occurrence of one or more peaks is not a serious problem in comparison with the 
experimental results, which predict a smooth plateau for vs A^:0in fact, the experimental 
data are relative to different samples with different composition, and any eventual peak 
would be smeared out on the average. The trend in vs N predicted by this RG analysis 
basically confirms the analogous result recovered via direct solution of the standard mean- 
field gap equation,EI even if Tc is never rigorously zero, except at the band edges. At variance 
with the experimental evidence, for uq = 0, has a steep increase at the band bottom, 
starting exactly at A^ = 0. We shall comment on this fact later on. 

In the case of such a symmetric interaction, = Uy, the exchange symmetry x ^ y, 
which is present in the hole dispersion relation ^(k) Eq. (|53|) and therefore in the Fermi 
surface, is restored even for the interaction. Thus, the kernel V{6,6') is now symmetric 
under the exchange 0^ ^ 6y, O'^ ^ d'y The eigenfunctions of the integral operator V are 
even or odd with respect to such an invariance transformation, and, as discussed at the end 
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of Section III, the scale invariant gap equation (|^) admits even (s wave) or odd {d wave) 
solutions for the gap function. 

Moreover, if uq = 0, there are just two eigenfunctions in the expansion for the gap, 
Eq. (^71), and at T = the gap equation admits the decoupled solutions: 

A(^) = const $i(^); A(^) = const ^2(0), (T = 0), (69) 

where the proportionality constants are fixed by the scale Aq defined by Eq. (|49|) and 
explicitly given by Eq. (|50|) : 

Ao = 2A4exp(A,-^(l,2)). (70) 

The eigenfunctions $1, $2 follow by direct solution of the system Eq. (0), which yields 

7y(l) = 7x(l), 72/(2) = -7a:(2), whence: 

^i{9)= const [cos9x + cos9y] s wave, 

$2(6')= const [cosOx ~ cosOy] c? wave. (71) 

The latter equations, together with Eq. ([7D|), fix the two symmetric solutions for the gap 
function. While in general a broken-symmetry solution may occur {e.g., a mixed s-d wavM), 
a pure s or d wave solution is always expected close to the transition points, which for T = 
are N = and ^ 0.5. In Fig. ^ we observe a cross-over around the peaked region of the 
eigenvalues: at the bottom of the band, /i = /ii, the most negative eigenvalue corresponds 
to the even eigenfunction $1, while for larger values of the most negative eigenvalue 
corresponds to the odd eigenfunction $2- The cross-over explains the occurrence of a double 
peak for Tc in Fig. and allows us to predict, at the transition point, the opening of an 
s-wave gap for < 0.2, and of a d-wave gap for A^ > 0.2. At T = 0, an intermediate 
broken-symmetry solution is awaited around the cross-over. All this is in agreement with a 
previous mean-field analysis.0 

Now let us switch on an on-site interaction, uq 7^ 0. In general, the cubic equation 
(|66|) may be solved analytically, but it is instructive to look a little closer at the case of a 
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symmetric inter-site interaction, = Uy. Again, the eigenfunctions of V must be even or 
odd with respect to the exchange x ^ y. 

The three hnearly independent functions Uq, Ux, Uy may be written in terms of the set 
{Uq, U±}, being U± = [Ux ± Uy]/\^. Since both Uq and f/+ are even, then the function f/_, 
which is odd, generates a one-dimensional invariant subspace for V, so that t/_ = $2 must 
be an eigenfunction for V, and its eigenvalue A (2) cannot be affected by the presence of an 
on-site coupling. Therefore, 

$2= const [cos9x — cos9x], (72) 
Xi2)= 2u4{U!) - {UxUy)]. (73) 

In fact, in the symmetric case, the cubic equation ( |66D factorizes as: 

(74) 

About the choice of the parameters Uq, u^, we may anticipate that the overall trend for 
Tc vs N is not relevantly affected by any reasonable change in their values, since it is a direct 
consequence of the chosen dispersion relation, as it is quite evident by comparison with the 
density of states pictured in Fig. ^. Nonetheless, the most interesting scenario shows up when 
there is a competition between on-site and inter-site interactions. Namely, for uq = —1.0, 
Ux = Uy = 1.0, the eigenvalues are reported in Fig. |^, and are characterized by some nice 
features: first of all, one only eigenvalue, A(0), is negative and leads the phase-transition; we 
may then notice that at the band bottom, — 0, the negative eigenvalue drops to a very 
small value at a finite filling A^ > 0; finally, the same eigenvalue goes to zero almost linearly 
inside the band for yU ~ 2.26, where a phase transition is expected even at T = 0. Actually, 
the eigenvalue never crosses the A = axis, as it is shown in the insert of Fig. ^. The 
eigenvalues A(0) and A(l), whose corresponding eigenfunctions share the same symmetry, 
cannot cross and repel each other, as predicted by a general theorem due to E. Wigner and 
J. von Neumann.0 However, we notice that a small negative eigenvalue is equivalent to a 
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mi) = ^{ui)+uxm) + {UxUy)]± 



^{u^)-uxm) + {u.Uy)) 



vanishing value since, due to the exponential in Eq. (p4D, the critical temperature becomes 
extremely low and thus negligible. 

Comparing with the case of a pure negative on-site interaction {ux = Uy = 0), when 
the eigenvalue coincides (up to a factor) with the density of states n of Fig. ^, we can 
assert that the presence of a non-zero nearest-neighbour repulsion reduces the range of the 
superconductive phase inside the band. This is not trivial since, for instance, if Ux < 0, the 
presence of an on-site repulsive interaction uq > does not produce any relevant effect on 
the phase diagram. In fact, reversing the sign of all the eigenvalues in Fig. and comparing 
with Fig. 1^, we notice that the presence of the d wave eigenvalue A(2), whose value is 
not affected by Uq, largely reduces the weight of an on-site repulsion against the inter-site 
attractive coupling. 

The progressive reducing of the superconducting phase with the increase of the inter-site 
repulsion is illustrated in Fig. |1^ and Fig. |ll[ Here, uq = —4.5 and Ux varies from 0.5 to 
3.5. Fig. [ly displays the negative eigenvalue while the critical temperature is reported in 
Fig. We notice that, as a consequence of the competition between uq and Ux, now at the 
bottom of the band the critical temperature drops to zero at a finite filling ^ 0, ranging 
from to 0.1, in qualitative agreement with the experimental data.0 

At the critical point the gap function is described by the even eigenfunction: 

^(0) =70 + 7a. (cos 0^. + cos^j,), (75) 

being (7o,7x,7y = 7x) the solution of the linear equations (|64D. 

Up to now we did not make any effort in order to justify the physical origin of the inter- 
action Eq. (|55D. This simple RG analysis seems to suggest the occurrence of a competition 
between an on-site negative Hubbard coupling and a repulsive nearest-neighbour interaction, 
on the basis of the comparison with the experimental data. We notice that, while several 
microscopic models predict the occurrence of a short-range attractive coupling,!! the nearest- 
neighbour repulsion could be justified by the presence of a long-range Coulomb interaction. 
Nonetheless, we must caution, since for a full analysis of the pairing interaction a larger 
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number of terms should be retained in its definition Eq. (^41), thus increasing the order of 
the couphng matrix, Eq. (pO|). 
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FIGURES 

FIG. 1. The 'tree' and the 'tadpole' diagrams are shown, which contribute to the interaction 
part of the action, Eq. at a tree level. 

FIG. 2. We show the three second-order one-loop diagrams, together with the factors arising 
from symmetry and statistics requirements. The only one which gives a relevant contribution is 
the BCS one.0 

FIG. 3. RG flow diagram. The flow lines for (/3A, A) have been numerically evaluated for 
different values of the unrenormalized eigenvalue Xq. The inverse temperature P is accordingly 
seen to renormalize either to zero or to some finite value, as soon as Xt is able to diverge at a finite 
t = i. The dashed line joins the critical points for the different values of A, and tends asymptotically 
to the axis /3 = as t — > oo. 

FIG. 4. The Fermi surfaces corresponding to a hole dispersion relation Eq. (^) are shown, 
in correspondence to increasing values of the chemical potential, /x = —2.1, — 1.8, 0.0, 3.0, 5.8. 
The box selects only the positive octant of the first Brillouin zone (0 < kx,ky,kz < tt), being 
F = (0, 0, 0), M = (vr, 0, 0), X = (vr, vr, 0), F = (vr, vr, vr). 

FIG. 5. The density of states n(;u), Eq. (|57D, is shown for the normal (non-interacting) Fermi 
system, as function of the chemical potential /x, /ii < /i < /is. It displays a pronounced, yet finite, 
maximum for /U = /i2 = —2.0. It is a token of the quasi-bidimensionality of the system, due to its 
layered structure, and would have been a real van Hove singularity if no inter-plane hopping had 
been considered (C = 0). 

FIG. 6. The total fraction of occupied states N{iJ,), Eq. (|58|), is shown for the normal 
(non-interacting) Fermi system, as functions of the chemical potential /i, /^i < /i < /U3. The 
normalization condition at the top of the band N^i^) = 1 is clearly fulfilled. 

FIG. 7. The two non-zero eigenvalues A(l) (dashed line) and A(2) (full line) are here displayed 
over /X, /Ui < /i < /i3, for the values of the parameters uq = 0.0, Ux = Uy = —1.0. 
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FIG. 8. We show here the ratio Tc/A for the values of the parameters uq = 0.0, Ux = Uy = —1.0, 
as a function of the fraction of occupied states N. 

FIG. 9. We show the eigenvalues A(0) (negative), A(l) (positive, dashed) and A(2) (full line) 
over Hi < < M3, in the symmetrical case corresponding to the values of the parameters 
no = —1.0, Ux = Uy = 1.0. The detail focuses on the region around fi ~ 2.262, where the two 
eigenvalues corresponding to the same symmetry (A(0) and A(l)) run closely without crossing. The 
eigenvalue A(0) < is seen to generate an instability. 

FIG. 10. We show the only eigenvalue giving rise to an instability, A(0), as a function 
of /X, /xi < /X < /X3, for the values of the parameters (from bottom to top): uq = —4.5, 
Ux = Uy = 0.5, 1.5, 2.5, 3.5. 

FIG. 11. We show the ratio Tc/A as a function of the fraction of occupied states iV(/x), for the 
values of the parameters (from top to bottom): uq = —4.5, Ux = Uy = 0.5, 1.5, 2.5, 3.5. 
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